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ABSTRACT 

Current constraints on dark matter density profiles from weak lensing are typically limited to 
radial scales greater than 50-100 kpc. In this paper, we explore the possibility of probing the 
very inner regions of galaxy/halo density profiles by measuring stacked weak lensing on scales 
of only a few tens of kpc. Our forecasts focus on scales smaller than the “equality radius” 
(R eq ) where the stellar component and the dark matter component contribute equally to the 
lensing signal. We compute the evolution of R eq as a function of lens stellar mass and redshift 
and show that R eq = 7-34 kpc for galaxies with M* = 10 9 - 5 -10 11 ' 5 M g . Unbiased shear 
measurements will be challenging on these scales. We introduce a simple metric to quantify 
how many source galaxies overlap with their neighbours and for which shear measurements 
will be challenging. Rejecting source galaxies with close-by companions results in a ~ 20 
per cent decrease in the overall source density. Despite this decrease, we show that Euclid 
and WFIRST will be able to constrain galaxy/halo density profiles at R eq with S/N > 20 for 
M* > 10 1o Mq. Weak lensing measurements at R eq , in combination with stellar kinematics 
on smaller scales, will be a powerful means by which to constrain both the inner slope of the 
dark matter density profile as well as the mass and redshift dependence of the stellar initial 
mass function. 

Key words: cosmology: observations - cosmology: large-scale structure of Universe - grav¬ 
itational lensing: weak 


1 INTRODUCTION 

According to the current cosmological framework, structure for¬ 
mation in the Universe is driven by the dynamics of cold dark 
matter. The collisionless gravitational collapse of dark matter over¬ 
densities and their subsequent virialization leads to the formation 
of dark matter halos with different masses and sizes. A variety of 
large scale cold dark matter numerical simulations (e.g., Dubinski 
& Carlberg 1991; Navarro et al. 1996b) have shown that the mass 
density distribution within halos is generally well described by a 
one parameter family, commonly referred to as the Navarro, Frenk, 
& White density profile (Navarro et al. 1997, hereafter NFW pro¬ 
file). 

The power-law slope of the dark matter distribution on scales 
of a few kpc to a few tens of kpc can provide clues to the nature of 
dark matter (Flores & Primack 1994; Moore 1994; Weinberg et al. 
2013). In particular, models with significant warm dark matter (e.g.. 


Maccio et al. 2012) or large self-interaction cross-sections (Spergel 
& Steinhardt 2000; Peter et al. 2013; Zavala et al. 2013) can result 
in a shallower density distribution in the core. However, baryonic 
physics also can significantly alter the distribution of dark matter on 
small scales, either by feedback (Navarro et al. 1996a; Mashchenko 
et al. 2006; Zolotov et al. 2012; Arraki et al. 2014) or simply by 
gravitational effects (adiabatic contraction: Blumenthal et al. 1986; 
Gnedin et al. 2004; Sellwood & McGaugh 2005). 

On radial scales below about one effective radius, the total 
mass profiles of galaxies transition from a dark matter dominated 
regime to a star-dominated regime. In addition, gas may represent 
a significant contribution in low mass galaxies. Disentangling the 
dark matter component from the stellar component on these scales 
is challenging. There are significant systematic uncertainties in the 
determination of galaxy stellar masses from the integrated light 
coming out of stars. Variations in the stellar initial mass function 
(hereafter IMF) and the low mass cut-off for star formation re- 
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suit in a factor of two uncertainty in stellar mass estimates (e.g., 
Barnabe et al. 2013; Courteau et al. 2014). Moreover, the IMF may 
vary with galaxy type and cosmic time (van Dokkum 2008; Con¬ 
roy et al. 2009; van Dokkum & Conroy 2010; Dutton et al. 2011a; 
Conroy & van Dokkum 2012; Smith et al. 2012). 

The dynamics of stars in the inner regions of galaxies (i.e. 
stellar kinematics) can be used to infer the total density profiles of 
galaxies on scales of about one effective radius. Using integral-field 
spectroscopy, the ATLAS3D collaboration obtained resolved stellar 
velocity dispersions of nearby elliptical galaxies to constrain stellar 
mass-to-light ratios (Cappellari et al. 2012). Another complemen¬ 
tary method is strong lensing, which provides fairly model indepen¬ 
dent constraints on the total mass within the Einstein radius. The 
combination of stellar kinematics and strong lensing has proved to 
be a very powerful approach to constrain the total density profile of 
galaxies on scales of about 3-9 kpc (e.g.. Sand et al. 2004; Koop- 
mans et al. 2006; Gavazzi et al. 2007; Jiang & Kochanek 2007; 
Auger et al. 2010a,b; Lagattuta et al. 2010; Dutton et al. 2011b; 
Newman et al. 2013a,b; Oguri et al. 2014; Sonnenfeld et al. 2014). 
Group scale strong lenses (with larger image separations) have also 
been used to probe the total density profiles of galaxies on ~ 10 - 
20 kpc scales (Kochanek & White 2001; Oguri 2006; More et al. 
2012). However, strong lensing systems are rare (at most a handful 
per square degree) and are primarily limited to massive early type 
galaxies at intermediate redshifts. In this paper, we explore the pos¬ 
sibility of using weak lensing measurements on small scales to po¬ 
tentially overcome these limitations and to probe the total density 
profiles of galaxies over a wide range in redshift and stellar mass 
(see e.g., Miyatake et al. 2013; More et al. 2014, for limits on stellar 
masses from weak lensing around BOSS galaxies). 

Our goal in this paper is to investigate how well future weak 
lensing surveys will be able to measure total density profiles at very 
small radial scales. In particular, we will focus on the transition 
point where the lensing signal is sensitive to both the stellar compo¬ 
nent and the dark matter component of the density profile. We call 
this transition scale the “equality radius”, noted hereafter as _R eq . 
Weak lensing measurements at the equality radius, in combina¬ 
tion with stellar kinematics on smaller scales, would be a powerful 
means by which to constrain both the inner slope of the dark matter 
density profile as well as the mass and redshift dependence of the 
IMF. To date, there have been few weak lensing measurements at 
r < R eq . Gavazzi et al. (2007) used Hubble Space Telescope imag¬ 
ing of 22 massive galaxies at z = 0.1-0.4 to measure three data 
points at r < R eq . The small number of background galaxies at 
these separations, however, results in large errors. Here we present 
predictions for how these types of measurements will improve as a 
function of galaxy mass and redshift for space-based weak lensing 
surveys such as Euclid (Laureijs et al. 2011) and WFIRST (Spergel 
et al. 2013). 

The measurement of a weak lensing signal at the equality ra¬ 
dius is inherently difficult. The equality radius is typically about 
20 kpc (See Section 3.1). This leaves only a tiny radial window in 
which to find and measure the shapes of background galaxies. Un¬ 
biased measurements of galaxy shapes in crowded environments 
poses yet another challenge. The light from the main lens (or light 
from neighbouring galaxies that correlate with the lens) may bias 
the shape measurements of source galaxies. The investigation of 
these biases is a topic of ongoing research, but will not be addressed 
in this paper. Instead, here we focus on more simple first-order 
questions. We will first investigate how the equality radius varies 
as a function of galaxy mass and redshift. Then, we will estimate 
how many source galaxies typically lie within this radial scale. We 


will also quantify how many source galaxies at R eq have very close 
companions and for which shape measurements will be challeng¬ 
ing. Finally, we would also like to understand whether there is a 
specific redshift and/or stellar mass range which is most suited for 
such measurements. 

This paper is organized as follows. After briefly presenting 
our theoretical framework in Section 2 , we investigate how R eq 
varies with stellar mass and redshift and explore various aspects 
that determine the signal-to-noise ratio (S/N) of weak lensing mea¬ 
surements at R eq in Section 3. The data analysed in our work is 
presented in Section 4. Our methodology is presented in Section 5. 
In section 5.2 we explore how cuts based on proximity effects im¬ 
pact the overall source density - these estimates are of general in¬ 
terest for all weak lensing studies including efforts to measure cos¬ 
mic shear. Finally, Section 6 shows our predictions for COSMOS 
(Scoville et al. 2007 ), WFIRST, and Euclid. We summarize our re¬ 
sults and present our conclusions in Section 7 . 

All radial scales are expressed in physical units. The projected 
transverse distance from a lens is denoted r whereas R^ rep¬ 
resents a three-dimensional distance. We define dark matter ha¬ 
los to enclose a spherical overdensity in which the mean den¬ 
sity is 200 times the mean background matter density, M 200b = 
| 7 rT? 2 oob 200 p b , where p b is the mean background density and 
J? 200 b is the halo boundary. For consistency with Leauthaud et al. 
( 2011 , 2012 ) we assume a WMAP5 (Hinshaw et al. 2009 ) cosmol¬ 
ogy with parameters = 0 . 258 , 17 a = 0 . 742 , Q. b h 2 = 0 . 02273 , 
n B = 0 . 963 , (j 8 = 0 . 796 , Ho = 72 kms _ 1 Mpc _1 . 


2 THEORETICAL FRAMEWORK 

In this section, we review the necessary theoretical background and 
introduce the concept of the “equality radius” which will be used 
extensively throughout this paper. 


2.1 Density Profiles of Stars and Dark Matter 

Weak gravitational lensing is the deflection of the path of light from 
distant galaxies due to the presence of mass over-densities along the 
line-of-sight. Weak lensing leads to both a distortion in the shapes 
of background galaxies and a magnification of galaxy fluxes and 
sizes. These effects are characterized respectively by the shear 7 
and the convergence k. In this work, we focus specifically on mea¬ 
surements of galaxy-galaxy lensing determined from the shear 7 . 

In the weak lensing regime, the average tangential ellipticity 
of background galaxies is related to the tangential component of 
the shear 7 t, which in turn is related to the excess surface density 
of the intervening mass 


7>w _ HrIAA) _ AM, (1) 

^crit ^crit 

where E(r) is the mean projected mass density within radius r and 
E(r) is the projected mass density at radius r. E cr i t is a critical 
density defined as: 


E 


crit 


c 2 Dos 
47tG D 0 iDi a ’ 


( 2 ) 


where D OB , D 0 \ and D\ s are respectively the angular diameter dis¬ 
tances from the observer to the source, from the observer to the 
lens, and from the lens to the source. 
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The excess surface density AS is additive and can be decom¬ 
posed into three terms representing contributions from dark matter, 
gas, and stars within galaxies: 


AE(r) = ASdm(r) + AS gaa (r) + AE ate ii a r(r) . (3) 


The gas component, AE gas , may have contributions from dif¬ 
ferent gas phases such as the cold interstellar medium (ISM) or the 
hot X-ray halo gas. As will be discussed further in Section 7, for 
most galaxies in our mass range, AE gaa is sub-dominant compared 
to the stellar and the dark matter contributions. For simplicity, we 
will neglect the AE gaa component in this work. 

To model the dark matter component, AEdm, we will assume 
the standard NFW profile: 


p{R 3 d) 


_^cPb_ 

(. R 3 d/R s ) (1 + R 3 d/R s ) 2 ’ 


(4) 


where R s is the characteristic scale radius of the halo. The charac¬ 
teristic density 5 C is given by 


g _ 200 _ cloob _ ^ 

3 lll(l + C 200 b) — c 200b/(l + C 200b) 

where the concentration parameter c 2 oo b is equal to the ratio 
Riaob/Rs- To compute c 2 oob as a function of M 20 ob and z, we 
adopt the concentration mass relation from Maccio et al. (2008). 
The AS profile of an NFW halo is given by 

AEdm(r) = r s 5 c p b g(r/R a ) , (6) 

where the function g(x) is a dimensionless profile (see Wright & 
Brainerd 2000). 

We model the stellar component using a Flernquist profile 
(Flernquist 1990). Flowever, over most of the radial range for which 
we can measure the shapes of background galaxies, the stellar com¬ 
ponent can be approximated by a simple point source term. For 
convenience we will adopt this approximation for most of our cal¬ 
culations. Under this approximation, AE ate ii a r is simply: 


AS a tell a r (t*) ~ \ . (7) 

nr Z 

Fig. 1 shows AS s teii a r computed using both the Flernquist 
profile and the point source approximation. 

Our goal in this paper is to study the radial scale at which the 
weak lensing signal becomes sensitive to the stellar mass of the 
lens sample. We define the “equality radius” as the radius at which 
AEdm(r) = AE a teii a r(r). This characteristic radius will be de¬ 
noted as R eq . On scales smaller than the equality radius, the weak 
galaxy-galaxy lensing signal will be dominated by the AE ate ii a r 
term and will hence directly probe M„. Fig. 1 gives an example of 
a typical AS profile and the equality radius R eq . 


2.2 Stellar-to-Halo Mass Relation 

To estimate M 200b from M* for each lens galaxy, we adopt 
the stellar-to-halo mass relation (SF1MR) from Leauthaud et al. 
(2012) (hereafter L12). In their study, L12 assume a lognor¬ 
mal distribution for the central galaxy conditional stellar mass 
function P(M* |A7 200b ). For central galaxies that reside in ha¬ 
los of mass M 200b , the average of their logarithmic stellar mass, 
/sHMR,(M 200 b), satisfies the following equation: 


lo SlO (/sHMr(^»)) — l°SlO (M 20 0b) 


= l°gio (Ml)+ /3login 



| (M./M., 0 ) s 

1 + (M*/M*,o) 1 


1 

2 ' 

(8) 


This functional form with parameters Mi, A/*,o, /?, S, and 7 is 
motivated from Behroozi et al. (2010). We refer the reader to Leau¬ 
thaud et al. (2011) and Behroozi et al. (2010) for the definitions of 
these parameters, notations, and how these parameters control the 
shape of SHMR. 

LI2 fit this model to the abundances of galaxies, their cluster¬ 
ing, and the weak lensing signal from the COSMOS field in three 
redshift bins (z 6 [0.22,0.48], [0.48,0.74], [0.74,1.00]). Their 
results suggest a very weak evolution of the global SFIMR from 
z = 0.2 to z = 1 , especially for low stellar mass galaxies with 
M* <2x 10 10 M 0 . We use the SFIMR parameters from L12 in 
each redshift bin. 

Note that this SFIMR is only valid for central galaxies. Ac¬ 
cording to the modelling results from LI 2 , about 20 per cent of 
galaxies in the COSMOS field are expected to be satellites. For our 
predictions, however, we neglect satellite galaxies and make the 
simplifying assumption that all lens galaxies are centrals. 

The average halo mass of galaxies of a given stellar mass 
can be computed using P(M 200 b |M*), which is related to 
P(M„ |A7 200b ) via Bayes’ theorem 


-P(M 200 b |M*) — 


oc 


P(M, | A/ 200 b )P(M 200b ) 
P(M„) 

P(M, |A/ 20 ob)n(M 20 ob) 1 


(9) 


where n(M 200b ) represents the halo mass function. For this anal¬ 
ysis, we will use the Tinker et al. (2008) halo mass function. We 
calculate the average halo mass of galaxies within a given stellar 
mass bin [Mi, M 2 ] as: 

Imi dM* f Q M 200 b P(M 20 0 b[M*)dM 200b 
200 b m,M 2 ] - — J .00 P(M 200 b |M*)dM 200b ■ 

( 10 ) 

The variation of A/ 200b as a function of lens stellar mass and red¬ 
shift is shown in Fig. 2. 


3 SMALL SCALE WEAK LENSING: MAIN EFFECTS 

Our goal is to predict the expected S/N of weak lensing measure¬ 
ments within the equality radius R eq . A variety of different effects 
will impact the expected weak lensing signal at this radial scale. 
We list these effects below. 

(i) -Req depends on the SHMR, the concentration-mass relation, 
and their evolution with redshift. 

(ii) The angular diameter distance and its dependence on red¬ 
shift determines the apparent area on the sky covered by a given 
R oq . For a fixed number of lens galaxies, a larger apparent area on 
the sky will correspond to a larger number of source galaxies and 
hence a higher S/N. 

(iii) For a fixed survey area, the number of lens galaxies will 
increase with redshift because the survey covers a larger comoving 
volume. 

(iv) The strength of the lensing signal (i.e, the shear) depends on 
the value of E cr i t which depends on the redshift of the lens and the 
source. For a fixed source redshift, z s , the strength of the lensing 
signal peaks around z s /2. 
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Figure 1 . Predicted AS profile for galaxies with log 10 (M*) ~ 11.3 and at 2 ~ 0.68. The blue dotted line represents ASj m assuming an NFW profile. 
The red dashed line represents AS ste iiar computed using a Hernquist mass model. The thin grey line shows AS Bte n ar computed using the point source 
approximation (Equation 7). The black solid line denotes AS which is the sum of AE f | ln and AS at eiiar- The green thin vertical line labelled “Rg o” shows 
the mean radius encompassing 90 per cent of the flux from lens galaxies. Another green thin vertical line labelled “R50” shows the half light radius of lens 
galaxies. The green thick vertical line shows the equality radius i? e q. The purpose of this paper is to predict the S/N of weak lensing measurements on radial 
scales smaller than R eq . 


(v) The strength of the shear around lens galaxies will be large 
at R cq . We need to be mindful of the regime where lensing is no 
longer weak. 

(vi) Magnification bias may alter the observed source density. 

(vii) The amount of “real estate” between the radius at which 
the light from a lens becomes insignificant and R eq depends on the 
sizes of lens galaxies. 

(viii) Unbiased measurements of shapes are difficult for galaxies 
that are blended or that are close to other galaxies so proximity 
effects are important. 


We compute R cq by equating the right hand sides of Equa¬ 
tions 6 and 7. Fig. 2 shows the variation of R eq as a function of 
lens stellar mass and redshift. At a given redshift, R eq increases 
with lens stellar mass. The increase in R eq with stellar mass is in 
contrast to the A/„ /Af 200b — Af 200b relation, which shows a dis¬ 
tinct peak at a preferred halo mass scale of M 2 oob = 10 11,8 Mq 
(e.g., Behroozi et al. 2010; L12; Rodrfguez-Puebla et al. 2013). 

This trend can be easily understood in the following manner. 
In the inner regions of halos (r < r B ), the profile function in Equa¬ 
tion 6, g(r/r B ) « 1. Therefore, we obtain: 


Some of these trends ((i)-(vi)) can be computed analytically. 
However (vii) and (viii) are non trivial. In particular, quantifying 
the effects of close pairs is especially difficult. For this reason, our 
predictions will be largely based on numbers drawn directly from a 
COSMOS weak lensing catalogue (see Section 4). Below, we dis¬ 
cuss each of these effects in further detail with an emphasis on de¬ 
veloping an intuition for how each effect drives our predicted S/N. 
We also clarify which components of our model are computed an¬ 
alytically and which components are drawn from real data. 


Rc q cc 


( M * 1 ) pJ 

V (M 2 00b) 3 ) 


( 11 ) 


For a fixed redshift, R cq is proportional to the first term in the above 
equation. We compute (M 200b ) from Equation 10, and find that the 
power law slope of the (Af 200b ) — AT* relation is approximately 
unity at a stellar mass of about log 10 (A7*) = 10.5, which implies 
the approximate scaling log 10 (f?. eq ) oc 0.34 x log 10 (Af*). This 
scaling is shown in Fig. 2 by a solid black line. It agrees with the 
increase of R cq with stellar mass that we observe. 


3.1 Equality radius 

As defined in Section 2.1, R eq corresponds to the radius where the 
contributions of the dark matter component and the stellar com¬ 
ponent to AS are equal. The S/N of weak lensing measurement 
on such small scales depends on the number of background source 
galaxies that lie within a projected distance r < R eq . Hence, a 
larger value of R eq should result in a larger S/N. 


3.2 Effective Area within R Bq 

We will use the Kron ellipse (which roughly encompasses 90 per 
cent of the light) to define the spatial extent of lens galaxies (see 
Section 5.2). We define the effective area, A e s, as the angular area 
between the outskirts of lens galaxies (as traced by the Kron el¬ 
lipse) and R eq - this is a measure of how much “real estate” is 
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log 10 (M*[M S ]) log lo (M.[M 0 ]) 



log lo (M,[M 0 ]) 


Figure 2. Left: variation of M 200b as a function of stellar mass and redshift based on SHMR. Coloured lines represent three different redshift bins for lens 
galaxies. Solid lines show ( whereas dashed line show (log 1 Q(M*)}. Differences between the solid and dashed lines are due to scatter in the SHMR 

relation. The green solid line corresponding to the redshift range (0.62 < z < 0.72) has been shifted horizontally by 0.03 dex for presentation purposes. 
Middle: evolution of R eq as a function of stellar mass and redshift. Coloured lines represent three different redshift bins for lens galaxies. Grey lines indicate 
the geometric mean of the semi-major and the semi-minor length of the Kron ellipses of lens galaxies. As a reference, the typical half light radius of lens 
galaxies at z ~ 0.4 is R e = 2.4 kpc for M* Ri 5.5 X 10 9 Mq and R e = 3.5 kpc for M* « 1.3 X 10 11 Mg. Right: evolution of the effective area within 
R aq as a function of stellar mass and redshift. 


available for weak lensing measurements at R cq . Fig. 2 shows how 
A e ff varies as a function of lens stellar mass and redshift. The red¬ 
shift dependence of A e g is driven by the angular diameter distance. 
A fixed value of R eq corresponds to a larger effective area at lower 
redshifts. Fig. 2 shows that massive galaxies at low redshift have 
the largest effective area and by consequence, the largest number 
of source galaxies per lens within the equality radius. 


In this paper, we are primarily interested in predictions for the 
S/N of weak lensing measurements for future weak lensing sur¬ 
veys. Flowever, we note that Equation 13 implies that another strat¬ 
egy for maximizing the S/N of small scale weak lensing measure¬ 
ments is targeted observations of very low redshift massive galaxies 
(Gavazzi et al. 2007; Okabe et al. in prep.) 


3.3 Geometric Effects due to E cr i t 

The strength of the lensing signal depends on the value of E cr ;t, 
which depends on the redshifts of both lens and source galaxies. 
To gain an intuition for how E cr i t drives S/N, let us consider a 
simple scenario in which we have a fixed number of lens galaxies 
N\ and a fixed high redshift source plane z s . The projected number 
density of sources is denoted as n s . The error on AE within a fixed 
physical aperture f? ap scales as cr AS oc E crit where N pails 

is the number of lens-source pairs. The number of lens-source pairs 
is simply: 


A/pairs = Ni x n s x nRl p /D 2 ol . (12) 

Putting these two together, the D a \ term cancels and the error on 
AE within a fixed physical aperture f? ap scales as: 


Dos 

D \ s 


(13) 


For a fixed source redshift, a lower lens redshift results in a larger 
D\ s , which corresponds to smaller error on AE. Therefore, from a 
pure geometric point of view, the S/N for weak lensing measure¬ 
ments within a fixed aperture will increase at lower redshifts. This 
may seem in contrast to the common intuition that the S/N of weak 
lensing measurements peaks at about a redshift of z s /2. The reason 
for this difference is simply that Equation 13 assumes a fixed num¬ 
ber of lens galaxies N \. In a weak lensing survey however, the num¬ 
ber of lens galaxies per unit comoving line-of-sight distance will 
decrease towards low redshifts (causing 7V pa i rs to also decrease). 


3.4 Amplitude of Shear on Small Scales 

As we try to push towards small radial scales, the weak shear as¬ 
sumption, |t|<< 1, may no longer be valid. Most of the shear 
measurement pipelines are designed with cosmic shear studies in 
mind and are not necessarily well tested for large shear values. For 
example, the GRavitational lEnsing Accuracy Testing 3 (GREAT3) 
competition (Mandelbaum et al. 2014) only tested 7 up to values of 
0.05. In this section we calculate the amplitude of the shear signal 
at the equality radius to determine whether or not this still corre¬ 
sponds to a weak lensing regime. 

At r = R cq , AE is simply: 

AE(i? eq ) = 2 x —±- . (14) 

The left hand panel of Fig. 3 shows the dependence of AE(f? eq ) on 
lens stellar mass and redshift. As discussed in Section 3.1, at fixed 
redshift, R eq increases with stellar mass but with a shallow slope of 
0.34. Therefore, AE(7? eq ) increases with stellar mass with a power 
law slope of 0.32. Similarly, at a fixed stellar mass, R eq decreases 
with redshift, leading to an increase in AE(f?. oq ) with redshift. 

We now estimate the typical value of 7 at r = R eq . For sim¬ 
plicity, we assume that all source galaxies are located at a fixed 
redshift z s = 2. The results are shown in the right hand panel of 
Fig. 3. The difference between the contour shapes in the left and the 
right panel are due to the dependence of E cr it on the lens redshift; 
for fixed z s = 2, E cr i t reaches a minimum at intermediate red¬ 
shifts. We find that galaxies with M* > 10 11 Mg, on average, will 
have shear values of order 0.05-0.01 at R cq . Hence, these types of 
measurements will require a careful calibration of shear measure¬ 
ments up to values of about 0.1. Biases in the shear measurement at 
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such high values are currently being investigated by the ARCLETS 
collaboration 1 . 

Moreover, in this shear regime, the shapes of source galaxies 
are described by a combination of the shear and the convergence, 
called the reduced shear (Schneider & Seitz 1995). The interpreta¬ 
tion of the small-scale weak lensing will require modelling of the 
reduced shear rather than the shear (that we investigate here). How¬ 
ever, the difference between the two is expected to be relatively mi¬ 
nor (5-10 per cent) on the radial scales and the halo mass scales 
considered in our work (e.g., Mandelbaum et al. 2006; Johnston 
et al. 2007; Leauthaud et al. 2010). 

An important point to note here is that for our calculations, we 
are considering the value of the shear for the average halo at fixed 
lens stellar mass. The type of measurements that we are proposing 
here are not possible at the hearts of rich groups or clusters of galax¬ 
ies at intermediate redshifts where the weak lensing approximation 
is clearly not valid. For massive galaxies, however, there is a large 
scatter in M 2 o 0b at fixed stellar mass (see e.g., More et al. 2011; 
L12). For example, at a fixed stellar mass of log 10 (M„) = 11, 
there is a 0.46 dex spread (a factor of 2.9) in halo mass, and a 0.7 
dex spread (a factor of 5.0) at log 10 (M*) = 10.5. Our calculations 
are valid for the average halo at fixed stellar mass - not for galaxies 
in rich groups or clusters. 

3.5 Proximity Effects and Blends 

To measure weak lensing signals, accurate and unbiased measure¬ 
ments of both shapes and photometric redshifts are critical. Such 
measurements are difficult when the angular separation between the 
lens and the source galaxy is small. In particular, the equality radius 
is merely a factor of 2-3 larger than the typical size of lens galax¬ 
ies (as traced by the Kron ellipse, see Fig. 2). If a source galaxy is 
located too close to a lens, then the light from the lens galaxy may 
contaminate the photometry as well as the shape measurement of 
the source galaxy. Biases in shape measurements will directly trans¬ 
late into biases in the measured weak lensing signal. Biases in the 
photometry contribute indirectly via biases in photometric redshifts 
which are required to compute E cr i t . The magnitude of these prox¬ 
imity effects depend on the size distribution of lens galaxies but 
also on the size distribution of neighbouring galaxies that correlate 
with the lens sample. We will use existing COSMOS ACS data to 
estimate the magnitude of these effects. We will devote large parts 
of Section 5 in order to quantify the magnitude of proximity effects. 


4 DATA 

As discussed in the previous section, among the various factors af¬ 
fecting the expected S/N of weak lensing measurements on small 
scales, the number of source galaxies lost due to proximity effects 
are the hardest to quantify analytically. It is primarily for this rea¬ 
son that we opt to use a real catalogue for the basis of our study 
instead of taking a semi-analytical approach such as the one ex¬ 
plored by Dawson et al. (2014). We will use the COSMOS ACS 
catalogue (described in the subsequent section) as our primary cat¬ 
alogue with which to conduct this analysis. The COSMOS field is 
small which means that our choice represents a trade-off between 
area and high resolution imaging but, as we will show further on in 

1 ARCLETS: http : //www . het . brown . edu/people/ian/ 
ClustersChallenge/ 


this paper, proximity effects are the largest determining factor for 
our predictions which motivates our choice. We now describe the 
data products that we use in greater detail. 

4.1 COSMOS Weak Lensing Catalogue 

The COSMOS program has imaged the largest contiguous area 
(1.64 degrees 2 ) with the Hubble Space Telescope (HST) using the 
Advanced Camera for Sun’eys (ACS) Wide Field Channel (WFC) 
(Scoville et al. 2007). The 5 a point-source limiting depth is Tab = 
27.2 mag in a O'.'24 diameter aperture (Koekemoer et al. 2007) and 
the size of point spread function (hereafter PSF) is 0.12" (Leau¬ 
thaud et al. 2007). This combination of depth and exquisite resolu¬ 
tion over a moderately wide area motivates our choice to use this 
as our primary data set. The details of the COSMOS weak lensing 
catalogue have already been described in ample detail elsewhere 
(Massey et al. 2007; Leauthaud et al. 2007; Rhodes et al. 2007; 
L12). The COSMOS weak lensing catalogue contains 3.9 x 10 5 
galaxies with accurate shape measurements, which represents a 
number density of 66 source galaxies per arcminute 2 . 

For shape measurements, the COSMOS weak lensing cata¬ 
logue uses the RRG algorithm (Rhodes et al. 2000). The exact de¬ 
tails of the shape measurement procedure are mostly un-important 
for this paper. We do however use information about when a shape 
measurement was possible with the expectation that the failure rate 
of galaxy shape measurement will increase for close lens-source 
pair configurations. The RRG COSMOS catalogue does not ex¬ 
plicitly flag source galaxies that have close companions. However, 
galaxies may fail to converge on a shape measurement if the RRG 
algorithm fails to determine a centroid (which typically occurs for 
blends). This is in contrast to the LENSFIT algorithm, for example, 
which explicitly rejects galaxies in close pairs whose isophotes are 
overlap at two-sigma level of the pixel noise (Miller et al. 2007, 
2013). 

To select galaxies with precise shape measurements, we adopt 
the same four selection cuts as Leauthaud et al. (2007). In addition, 
we restrict ourselves to galaxies that are detected in the COSMOS 
Subaru catalogue (See Section 5.3) and with 0 < o me as < 0.2, 
where <r meaB represents the shape measurement error from Leau¬ 
thaud et al. (2007). We also remove a small fraction of galaxies for 
which SourceExtractor (Bertin & Arnouts 1996) failed to measure 
the KRON_RADIUS (see Section 5.2 for the definition of this pa¬ 
rameter). These additional cuts remove 29 per cent of galaxies from 
the COSMOS weak lensing catalogue, leaving 2.7 x 10 5 galaxies 
(45 galaxies per arcminute 2 ) with shape measurements and which 
are detected by the COSMOS Subaru observations. 


4.2 COSMOS Photo-z Catalogue 

Photometric redshifts (hereafter, “photo-zs”) are necessary in order 
to separate foreground and background galaxies. Ideally, our photo- 
zs would be determined from multi-band imaging matched in both 
resolution and depth to the COSMOS F814W imaging. However, 
the COSMOS photo-zs are measured from ground-based imaging 
with a PSF that is larger by about a factor of 10 compared to the 
ACS data. In Section 5, we study how this photo-z matching affects 
our predictions. 

For this paper, we use the COSMOS photo-z catalogue ver¬ 
sion 1.8 presented in Ilbert et al. (2009). These photo-zs have been 
derived using a x 2 template fitting method using over 30 bands of 
multi-wavelength data from UV, visible near-IR to mid-IR, and also 
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Figure 3. Left: variation of AS at _R eq as a function of lens stellar mass and redshift assuming source galaxies at z a = 2. AS(f? eq ) monotonically increases 
with stellar mass and redshift. Right: same as the left panel but for the shear 7 . At fixed stellar mass, 7 reaches maximum values at 2 ~ 0.6. 



have been calibrated with large spectroscopic samples from VLT- 
VIMOS (Lilly et al. 2007) and Keck-DEIMOS. 

4.3 Stellar Mass Estimates 

Stellar mass estimates are required in order to predict AS. We use 
the same stellar mass estimates as L12 and refer the reader to L12 
and Bundy et al. (2010) for further details. 

Stellar mass estimates are based on PSF-matched 3”0 diam¬ 
eter aperture photometry from the ground-based COSMOS cata¬ 
logue (filters u*, Bj, Vj, g + , r + , i + , z + , K s ) (Capak et al. 2007; 
Ilbert et al. 2009; McCracken et al. 2010). The depth in all bands 
reaches at least 25th magnitude (AB) except for the 7C s -band, 
which is limited to K s < 24 [mag]. Stellar masses are derived us¬ 
ing the Bayesian code described in Bundy et al. (2006), which uses 
Bruzual & Chariot (2003) population synthesis code and assume 
a Chabrier IMF (Chabrier 2003) and a Chariot & Fall (2000) dust 
model. The stellar mass completeness is defined from magnitude 
limits K s < 24 [mag] and /smw < 25 [mag]. For the redshift 
range that we are interested in here (0.2 < 2 < 1.0), the stellar 
mass completeness ranges from M, = 10 816 Mq to 1O 9,98 M 0 . 


5 METHOD 

Here we describe the method that we use to compute the predicted 
S/N for AE(r < R eq ) for weak lensing data sets with COSMOS 
like quality. For a given set of foreground lens galaxies, the main 
ingredients necessary to make this prediction are the number of 
lens galaxies, the number of background galaxies with shape mea¬ 
surements as a function of distance, the shape noise for these back¬ 
ground galaxies, and the redshift distribution of source galaxies. 
We describe how we extract these quantities from the COSMOS 
ACS catalogue. 

5.1 Lens and Source Samples 

There is some arbitrariness involved in the definition of the fore¬ 
ground lens sample used to perform a S/N investigation. The pre¬ 
dicted value of the S/N depends on the number of galaxies in the 
lens sample, which depends on the choice of the binning scheme. 


For example, equal-sized redshift bins will provide a larger volume 
for higher redshift bins, which results in a larger number of galax¬ 
ies in the higher redshift bins. In our work, we opt to divide our 
galaxies into 7 equal comoving volume bins in redshift and 6 log¬ 
arithmically spaced stellar mass bins. In each redshift bin, we only 
consider bins that are complete in stellar mass (Section 4.3; L12). 
The characteristics of our lens samples are provided in Table 1. 

Despite our choice of redshift bins with equal comoving vol¬ 
umes, the number of lens galaxies within each stellar mass bin fluc¬ 
tuates due to sample variance. We would like to minimize the ef¬ 
fects of such variations on the calculated S/N. We make the assump¬ 
tion that the total stellar mass function of galaxies does not strongly 
evolve at 2 < 1 (motivated by Bundy et al. 2010; Ilbert et al. 2010; 
L12; Moustakas et al. 2013). We fit a constant to the number of 
galaxies in each stellar mass bin, assuming Poisson statistics. The 
value of these constants are provided in Table 1 as “mean number”. 
We use this mean number in order to generate our predictions. 

As our input catalogue for source galaxies, we use the COS¬ 
MOS ACS weak lensing catalogue. As described in Section 4.1, 
we restrict ourselves to source galaxies that pass a combination of 
three quality cuts as well as four quality cuts described in Leau- 
thaud et al. (2007). These galaxies form our source catalogue. 


5.2 Impact of Proximity Effects on the Overall Source 

Density 

In this section we explore how cuts based on proximity effects im¬ 
pact the overall source density. These estimates are of general inter¬ 
est for all weak lensing studies including efforts to measure cosmic 
shear for example. 

Galaxies in close pairs may have biased shear estimates, but 
the exact level of such bias or how it varies as a function of the dis¬ 
tance between close pairs, remains poorly understood. As a result, 
most current weak lensing pipelines do not necessarily have a well 
justified criterion to determine which galaxies have shapes that may 
be biased because of neighbouring galaxies. For example, the ACS 
lensing catalogue that we use in this paper makes no stringent cuts 
to remove source galaxies in close-pair configurations, even though 
these galaxies may have biased shear estimates. On the other hand. 
Miller et al. (2013) take a more conservative approach while con¬ 
structing the CFHTLS weak lensing catalogue (see Section 4.1). 
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Redshift 

[ 9.37,9.74 ] 

[9.74,10.11 ] 

log 10 M, [Mq] 

[10.11,10.49] [10.49.10.86] 

[10.86, 11.23] 

[ 11.23 , 11.60] 

[0.30,0.50] 

2477 

2000 

1690 

1404 

731 

162 

[ 0.50,0.62 ] 

1688 

1185 

1110 

878 

494 

113 

[0.62,0.72] 

X 

1476 

1209 

1144 

750 

212 

[0.72,0.80] 

X 

1669 

888 

767 

479 

137 

[0.80,0.87] 

X 

X 

1436 

1458 

943 

266 

[0.87,0.94] 

X 

X 

1107 

993 

622 

173 

[ 0.94 , 1.00 ] 

X 

X 

1185 

1058 

737 

241 

Mean number 

2007.8 

1526.4 

1188.6 

1049.9 

645.9 

171.9 


Table 1. Number of lens galaxies from our COSMOS catalogue in bins of stellar mass and redshift. Bins that are incomplete in terms of stellar mass are shown 
as X. The numbers in the bottom row are the mean values after the correction of sample variance described in Section 5.1. 


However, both of these choices remain subjective and do not study 
the amount of shear bias in close pairs. Quantitative investigation 
of shear bias is a topic of active research, but is beyond the scope 
of this paper. Instead, we use a simple and easily reproducible cri¬ 
terion for identifying close pairs and we explore how the predicted 
S/N for small-scale lensing measurements varies for source galaxy 
selections that are more or less conservative. 

Our scheme to identify close pairs of galaxies is based on 
the Kron parameters from SourceExtractor (Bertin & Arnouts 
1996) (hereafter SEXTRACTOR), which generalizes the luminosity- 
weighted radius of Kron (1980) to elliptical apertures. Our choice 
is motivated by the fact that Kron parameters are often available in 
imaging catalogue, which makes our criterion easily reproducible. 
We note, however, that our procedure will not identify blends at a 
fixed isophotal level. While an isophotal definition of blends would 
have been preferable, it was not possible to calculate such a crite¬ 
rion given the catalogues available to us. 

The Kron parameters consist of AJMAGE, BJMAGE, 
KRON JIADIUS, and THETA.WORLD. These parameters can 
be used to define a “Kron ellipse” for each galaxy that roughly 
encompasses 90-95 per cent of the flux. For SEXTRACTOR, 
KRON_RADIUS is computed by multiplying the first moment ra¬ 
dius by a factor of 2.5, then the major axis and the minor axis of 
the Kron ellipse are given by AJMAGE x KRON JIADIUS and 
BJMAGE x KRON_RADIUS respectively. The position angle 
of the Kron ellipse, THETA_WORLD, is measured between the 
right ascension and the direction of the semi-major axis. 

We place Kron ellipses around all galaxies and identify galax¬ 
ies with overlapping Kron ellipses. Our baseline predictions reject 
all source galaxies with overlapping Kron ellipses. In order to make 
this rejection criterion more or less conservative, we simply in¬ 
crease or decrease the size of the ellipse by multiplying the ma¬ 
jor and minor axis by a constant factor /. This procedure is illus¬ 
trated in Fig. 4. The subsample of galaxies rejected by our criterion, 
when the multiplicative factor is /, will be denoted by Rf. Those 
galaxies which are rejected when f = f 2 but not rejected when 
/ = /i < f 2 , will be denoted by We use / = 1-0 for our 

fiducial set of predictions and we explore how our predictions vary 
for different values of /. Table 2 shows the number and fraction of 
source galaxies that are rejected for various values of /. 

Table 2 shows that the overall source density is very sensitive 
to proximity cuts. Rejecting all source galaxies that have a Kron el¬ 
lipse that overlaps with a neighbouring galaxy leads to a 20 per cent 
decrease in the overall source density. Adopting more conservative 
values of / = 1.2 will lead to a 30 per cent decrease in the overall 
source density. Understanding how neighbouring galaxies impact 



Figure 4. Two neighbouring galaxies with Kron ellipses overlaid. For the 
blue ellipses, the Kron ellipse has been scaled by a factor of / = 1.2. For 
the red ellipses, the Kron ellipse has been scaled by a factor of / = 0.8. 
These galaxies are rejected from our source catalogue when / = 1.2 but 
are not rejected when / = 0.8. The bright compact object is a star and does 
not belong in the source catalogue. 


-R[0,0.5] 

^[0,0.8] 

R[o,l.o] 

R[0,1.2] 

■R[0,1.5] 

Number 16993 

43182 

61023 

80473 

111577 

Fraction(%) 6.3 

16.0 

22.6 

29.8 

41.3 


Table 2. The number and fraction of galaxies identified as close pairs in the 
source catalogue for different values of /. Our default scheme is / = 1. 
In this case, 22.6 per cent of galaxies are rejected from the overall source 
catalogue. 


shear bias and which source galaxies need to be rejected is clearly 
of importance for all weak lensing studies, not just the particular 
science application discussed in this paper. 

Fig. 5 shows an example of source galaxies selected by our 
method. In this figure, we show small cutout images around three 
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massive lens galaxies at z ~ 0.25. Red ellipses indicate galax¬ 
ies that have overlapping Rron ellipses. To predict the S/N of AE 
measurements within R eq , we only use source galaxies with green 
ellipses. 

5.3 Impact of Ground-Based Photometric Redshifts on 
Source Counts 

In addition to shape measurements, weak lensing also requires 
photometric redshifts. The assignment of photometric redshifts to 
galaxies in close-pair configurations is non-trivial. Photometric red¬ 
shifts are often derived from ground based imaging. For example, 
both Euclid (Laureijs et al. 2011) and WFIRST (Spergel et al. 2013) 
will need to complement their space-based imaging with ground 
based photometry. Because ground based imaging typically has a 
PSF of about one arcsecond, it may be difficult to derive accurate 
photometry for source galaxies at R cq . 

We use the COSMOS photo-z catalogue to quantify how many 
source galaxies we lose due to this additional photometric redshift 
requirement. The detection of sources in this catalogue was carried 
out on a combined CFHT i* and Subaru i + image (with the original 
PSF of 0'.'95) (Capak et al. 2007). Galaxy colours were measured 
with a 3V0 aperture after PSF homogenization. 

We now investigate the conditions under which shape and pho¬ 
tometric redshift measurements fail in the COSMOS catalogues. 
Let us consider a galaxy A with a close-by companion galaxy B. 
We would like to know when shape and photo-z measurements typ¬ 
ically fail for galaxy A as a function of the distance from B and as 
a function of the brightness ratio between galaxies A and B. To an¬ 
swer this question, we bin all of the galaxies in our ACS catalogue 
by magnitude and refer to it as our “primary” sample. For every 
galaxy in this sample, we identify all its neighbours in the entire 
ACS catalogue 2 . We use a criteria based on Kron ellipses in order 
to scale distances between the primary and its neighbours. Pairs 
of primary-neighbour galaxies whose Kron ellipses overlap when 
their major and minor axes are scaled by a factor fi < f < fz, are 
labelled R\f 1 j 2 1. If there is more than one overlapping neighbour, 
we use the brightest one (this occurs in about 6 to 10 per cent of 
cases depending upon the value of /). 

Fig. 6 shows how often shapes and photo-zs are measured for 
primary galaxies as a function of the distance and magnitude dif¬ 
ference with the brightest neighbour. The two columns correspond 
to different primary samples (the fainter primary sample is on the 
right). The top row corresponds to galaxies in close pairs where the 
separations are small (i?[o.5,o.8])- We investigate close-pair galaxies 
with successively larger separations in the middle (i?[o.8,1.0]) and 
the bottom row (fin. 0 , 1 . 3 ])- Fig. 6 shows that shape and photo-z 
measurements are more likely to fail for galaxies which have bright 
companions. Comparing the left and the right hand panels shows 
that this effect is more severe for faint galaxies. Comparing the dif¬ 
ferent rows it is clear that the effect becomes smaller as we consider 
galaxies in close pairs but with larger separations. 

In this test, we only consider when a galaxy has a photomet¬ 
ric redshift measurement, but we do not evaluate whether or not 
photo-zs for galaxies with close by companions have a larger bias 
or a larger fraction of catastrophic errors. In addition, here we also 
only consider when a shape measurements has been possible, not 
whether or these shape measurements are biased. Obviously, these 
are aspects that need to be evaluated more critically in future work. 

2 These neighbours could themselves be a part of the primary sample. 


Bearing these caveats in mind, Fig. 6 shows that shape measure¬ 
ments in the COSMOS catalogue are a more stringent requirement 
compared to photo-zs. 

5.4 Proximity Effects in the Number of Source Galaxies as a 
Function of Transverse Distance from Lenses 

We investigate how blending affects the number of lens-source 
pairs as a function of transverse separation, r. The upper panels 
of Fig. 7 show the cumulative number of observed lens-source 
pairs, (Vpairs(< r), in two different stellar mass bins and at one 
fixed redshift bin. Symbols with different colours correspond to all 
HST detected galaxies, those with photometric redshifts, and those 
that pass our source galaxy cuts. Note, however, that we have not 
yet imposed any proximity cuts on the source catalogue - galax¬ 
ies with overlapping Kron ellipses are still included at this stage. 
For all samples in this figure, we also impose a magnitude cut at 
i < 26. Note that we also do not yet impose a photo-z cut to sep¬ 
arate foreground and background objects - we simply consider all 
pairs along the line-of-sight. 

Fig. 7 shows that the overall number of lens-source pairs at 
large separations (~ 40 kpc) decreases as we consider higher mass 
lens galaxies (indicated by the magenta vertical thin line). How¬ 
ever, high mass lens galaxies have a larger number of lens-source 
pairs within R cq (indicated by the magenta horizontal thin line). 
This is consistent with our expectations from Section 3.1: both R eq 
and Aefiect are larger for more massive galaxies. Because massive 
galaxies have a larger effective area where source galaxies can be 
found, they are more suited for detecting the small scale weak lens¬ 
ing signal. 

The bottom panels of Fig. 7 show the ratio of /V pa i rs (< r) to 
the number of pairs expected from a simple power law extrapola¬ 
tion (with slope —2) of iV pa i rs (< r) from larger scales (thin grey 
line). These panels show that although a number of galaxies are 
detected at r < R eq , our photo-z and shape measurement require¬ 
ments bring these numbers down significantly. Also, as discussed in 
the previous section, Fig. 7 shows that we lose more source galax¬ 
ies due to the shape measurement requirement than the photo-z 
requirement. Note that all HST pairs in the higher mass lens bin 
outnumber the large-scale extrapolation on scales ~ 20kpc. This 
bump feature indicates that galaxies are clustered with lens sam¬ 
ples along the line-of-sight. These galaxies can be misinterpreted 
as source galaxies in our S/N predictions due to the photo-z errors. 
In Section 6.2, we will discuss such potential confusion in source 
selection due to correlated galaxies. 

In addition to these effects, galaxies will also be rejected from 
our source catalogue due to the proximity cuts defined in Sec¬ 
tion 5.2. We are interested to know if source galaxies are flagged 
as overlaps at R aq because they overlap with the lens galaxy un¬ 
der consideration, or because they overlap with non-lens galaxies.. 
Fig. 8 quantifies the relative importance of the two effects. For low 
stellar mass galaxies, we find that the majority of source galaxies 
lost within R cq are due to overlap with the lens galaxy itself. At 
higher stellar masses we find that only ~30 per cent of the source 
galaxies flagged as overlaps are colliding with the primary lens 
sample. The remaining ~70 per cent overlap with either nearby 
galaxies that correlate with the primary lens sample, or with galax¬ 
ies that are spatially uncorrelated with the lens sample. Again in 
this figure, we consider all pairs along the line-of-sight without the 
photo-z cut. 

Finally, to select only source galaxies located at higher red¬ 
shifts than lens galaxies, we impose a photo-z cut z s > 21,max 
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Figure 5. The three large galaxies at the centre of each postage stamp are examples of lens galaxies from the COSMOS survey at z ~ 0.25 and with 
log i q M*[Mq] ~ 11.1. Green ellipses show Kron ellipses (here we use f = 1). Blue circles show R eq . Galaxies with overlapping Kron ellipses are 
identified by red ellipses - these are rejected from our source catalogue. No circles or ellipses are drawn for stars. 


where egmax represents the upper redshift limit of the lens bin un¬ 
der consideration [zi, min, zi.max]- The combination of the photo-z 
cut, the shape measurement cut and the proximity cut results in 
the cumulative number distribution of lens-source pairs shown in 
Fig. 9. Symbols with different colours correspond to source galax¬ 
ies that pass only photo-z cut and shape measurement cut (in red), 
those that pass also proximity cut with f = 1 (in green), and those 
that pass also proximity cut with =1.3 (in blue). In the fiducial 
case (/ = 1), the fraction of source galaxies within R eq that pass 
the photo-z cut, but are rejected due to a combination of the shape 
measurement cut and the proximity cut is over 80 per cent. We will 
use the A r pa i ra (< R eq ) in order to estimate the error on the weak 
lensing measurement within R cq for every stellar mass and redshift 
bin in our sample. We compute the weak lensing signal at the aver¬ 
age separation between the lens-source pairs using Equations 6 and 
7. 


5.5 7Vp a i rs (< R eq ) : Number of Lens-Source Pairs Within the 
Equality Radius 

We compute the total number of source galaxies that pass the 
photo-z cut, the shape measurement cut, the proximity cut, and 
which have r < R oq for each of our lens bins. The values for 
Adairs(< -Rcq) for / = 1 are given in Table 3. As can be seen 
from Table 3, the overall number of lens-source pairs for any given 
bin is small, only of order ~ 100. The difference between Tables 1 
and 3 agrees with the results our conclusions from Section 5.4 - 
low-mass lens galaxies have almost no source galaxies within R eq 
due to the small size of R oq . 


5.6 Error on AS and Signal to Noise Ratio 

Following Leauthaud et al. (2007), the error on the shear for each 
source galaxy is estimated as a combination of intrinsic shape noise 
(<r int = 0.27) and shape measurement error: 


2.2 
^int ' 0" meas . 


(15) 


The optimal weight for each source galaxy is given by 


Wi 


(Sc 


(16) 


and the estimated error on AS in the COSMOS survey is given by 


COSMOS 


^pairs ^teq) 


E 


- 1/2 


(17) 


The sum runs over all lens-source pairs with r < R oq in a given 
stellar mass and redshift bin. 

To scale our error estimate to survey areas larger than COS¬ 
MOS, we assume that the distribution of source galaxy redshift and 
shape measurement errors in COSMOS are representative and we 
simply scale the estimated errors according to 


survey COSMOS 


AcoSMOS 


A s 


(18) 


where Acosmos = 1-64 deg 2 . Finally, combining AS and has 
(Equations 3 and 17), the predicted S/N in each stellar mass and 
redshift bin for COSMOS is simply 


S/N 


AE 


(19) 


5.7 Predictions for COSMOS 

Using Equations 17 and 19, we compute the expected S/N for one 
weak lensing data point at r < R eq in the COSMOS ACS cat¬ 
alogue. The results are shown in Fig. 10. We find that the S/N 
for this type of measurement is maximized for massive galaxies 
at low redshifts, in agreement with our intuitive expectation from 
Section 3.4. We also show the effects of making more or less con¬ 
servative source galaxy selections by varying our / factor. A larger 
value of / means that we reject more source galaxies in close-pair 
configurations. We find that varying / between 0.8 and 1.3 only has 
a relatively minor impact on the predicted S/N. This implies that a 
more conservative selection of sources galaxies does not necessar¬ 
ily have a strong impact on the S/N. 
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Figure 6. Number of galaxies in COSMOS as a function of the magnitude difference (may (primary galaxy) - maglneighbour galaxy)) between galaxies in 
close-pair configurations. Left panels: primary galaxies with magnitudes S [24.0, 24.5]. The distance to the closest neighbour increases from top to bottom 
(If'try. 0 . 8 ] i s a small separation and -R[i. 0 , 1 . 3 ] is a large separation). Blue histograms represents all galaxies in the COSMOS catalogue. Green histograms 
(“Matched to ground”) represent primary galaxies that have a photometric redshift measurement. Red histograms (“With shape”) represent primary galaxies 
that have a shape measurement and that pass our quality cuts as discussed in Section 4.1. Right panels: same as left panels but for primary galaxies with 
magnitudes 6 [25.0, 25.5]. This figure shows that both shape and photo-z measurements fail for source galaxies with bright companions. For example, panel 
d) shows that shape measurements fail for about 50 per cent of galaxies with magnitude £ [25.0, 25.5] and which have a close-by companion of equal 
magnitude (A mag = 0). This failure rate increases sharply for brighter companions (Amag > 0). Panel e) however shows that at rnagG [24.0, 24.5] and 
for 1.0 < / < 1.3, most galaxies have both shape and photo-z measurements. A key point to note from this figure is that the red histograms are typically 
lower than the green: shape measurements in the COSMOS catalogue are a more stringent requirement compared to photo-zs. 


6 PREDICTIONS FOR EUCLID AND WFIRST 

In this section, we present our predictions for the S/N for small 
scale weak lensing measurements from future space-based surveys 
Euclid and WFIRST. 


6.1 Next Generation Weak Lensing Space Based Surveys 

Euclid is a space mission under development for an expected 
launch in 2020 (Laureijs et al. 2011). Euclid consists of a 1.2 m 
Korsch telescope and two instruments, the VIS (visible imager) and 
NISP (Near Infrared Spectrometer and Photometer) and is designed 
to study the properties of dark energy and dark matter and to search 


for evidence of modified gravity by using weak gravitational lens¬ 
ing and galaxy clustering. Euclid will perform a 6 year imaging 
and spectroscopic survey over the lowest background 15000 deg 2 
of the extragalactic sky. Visible imaging in a single wide riz fil¬ 
ter will provide shapes for about 1.6 billion galaxies; near infrared 
imaging in Y, J, H bands combined with ground-based photometry 
will enable high precision photometric redshifts for source galax¬ 
ies. Euclid will have a ~0.2 arcsecond PSF and will measure galaxy 
shapes with a source density of ~30 galaxies per arcmin 2 . 

The Wide Field Infrared Survey Telescope (WFIRST) is a 
NASA mission under study for possible launch in ~ 2024 (Spergel 
et al. 2013). WFIRST consists of 2.4m telescope a Wide Field Im¬ 
ager with 18 4kx4k near infrared (NIR) detectors with 0.1 arc- 
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Radius [kpc] Radius [kpc] 


Figure 7. Upper panels: cumulative number of lens-source pairs, -/V pa i rs (< r), as a function of the transverse distance from lens galaxies at z £ [0.62, 0.72]. 
Left panels correspond to lens galaxies with log 10 M* E [10.11,10.49]. Right panels correspond to lens galaxies with [10.86,11.23]. Blue dashed vertical 
lines indicate the average semi-major and semi-minor axis of Kron ellipses for lens galaxies. The thick red vertical line denotes R eq . Green dots indicate all 
galaxies from the COSMOS weak lensing catalogue with i < 26. Red dots represent galaxies with shape measurements. Blue dots represent galaxies that 
are matched to the ground-based photo-z catalogue. Note that no photo-z cuts have been applied to limit source galaxies to background galaxies (these cuts 
will be applied in Fig. 9. Magenta vertical thin line shows the location of r = 40 kpc, and magenta horizontal thin lines respectively show iV pa i rs (< 40kcp) 
and AT pa i rs (< R eq ). Lower panels: Af pa j rs (< r) divided by the expected number based on an extrapolation from larger scales (iVextra)- The lower panels 
show that proximity effects starts to influence the source galaxy counts at r ~ 40 kpc (roughly 1 to 5 times larger than R eq ) and photo-z matching and shape 


iVn 


s(< r). 


Redshift 

[ 9.37,9.74 ] 

[9.74,10.11 ] 

logic M*[M Q ] 

[10.11,10.49] [10.49,10.86] 

[ 10.86, 11.23] 

[ 11.23 , 11.60] 

[0.30,0.50] 

1 

14 

66 

197 

178 

51 

[0.50,0.62] 

0 

1 

21 

70 

84 

29 

[0.62,0.72] 

X 

1 

21 

60 

87 

48 

[0.72,0.80] 

X 

1 

12 

37 

52 

24 

[0.80,0.87] 

X 

X 

17 

55 

91 

39 

[ 0.87,0.94 ] 

X 

X 

6 

31 

43 

27 

[0.94, 1.00] 

X 

X 

7 

40 

57 

25 


Table 3. The number of lens-source pairs for r < R eq for each of our lens bins. Source galaxies which have Kron ellipses that overlap with any nearby galaxy 
(/= 1) are rejected from the source catalogue before computing these numbers. The photo-z cut to seperate background and foreground source galaxies has 
also been applied before compute these numbers. 


© 0000 RAS, MNRAS 000, 000-000 









































Weak leasing on 20 kiloparsec scales 


13 




Radius [kpc] 


Radius [kpc] 


Figure 8. Similar to Fig. 7 but here we only show galaxies which have a Kron ellipse that overlaps with a neighbouring galaxy (/ = 1). Upper panels: 
cumulative number of galaxies identified as close pairs as a function of the transverse distance from lens galaxies. We distinguish overlaps that occur between 
a source and a lens from the sample under consideration from overlaps that occur between a source and some other non-lens galaxy. Magenta dots indicate 
the cumulative number of all source galaxies that are flagged are overlaps. Green dots represent source galaxies that have a Kron ellipse that overlaps with 
the Kron ellipse of a lens galaxy. The overall amplitude of A r pa ; rs (< r) on larger radial scales is dominated by overlap between source galaxies that are 
spatially uncorrelated with the lens sample. On the other hand, the radial dependence of JV pa i ra (< r) on smaller radial scales is mainly driven by overlap with 
the primary lens sample and with galaxies that are correlated with the lens galaxies. Note that no photo-z cuts have been applied to limit source galaxies to 
background galaxies (these cuts will be applied in Fig. 9). Lower panels: fraction of source galaxies flagged as overlaps and for which the overlap occurs with 
a lens galaxy from the sample under consideration. Panel c) shows that for the lower stellar mass bins, overlaps mainly occur at R eq because source galaxies 
tend to overlap with galaxies from the primary lens sample. Panel d) on the other hand, shows that for the higher stellar mass bins, only ~30 per cent of the 
source galaxies flagged as overlaps are colliding with the primary lens sample. The remaining ~70 per cent overlap with either galaxies that are spatially 
correlated with the lens sample or with galaxies that are uncorrelated with the lens galaxies. 


second pixels, an integral field unit (IFU) spectrograph, and an 
exoplanet coronagraph. WFIRST is designed to perform NIR sur¬ 
veys for a wide range of astrophysics goals, including weak lens- 
ing. WFIRST s weak lensing survey will cover 2400 deg 2 and make 
galaxy shape measurements in 3 NIR bands of about 500 million 
galaxies over the course of two years in WFIRST s primary mission 
of 6 years. WFIRST will have a PSF of 0.1 arcseconds (before pix- 
ellization) at 1 micron and about 0.2 arcseconds at 2 microns and 
will measure galaxies shapes with a source density of ~54 galaxies 
per arcmin 2 . 

6.2 Predictions for one bin at r < R eq 

To make our predictions for Euclid, we assume a source den¬ 
sity of 30 galaxies per arcmin 2 over 15000 deg 2 . To mimic this 
source density, we simply apply a detection signal-to-noise cut on 


our COSMOS catalogue. After applying a signal-to-noise cut that 
yields 30 galaxies per arcmin 2 , we rederive N^clid (< f or 
Euclid. 

WFIRST is expected to achieve a source density of ~54 galax¬ 
ies per arcmin 2 over 2400 deg 2 . The COSMOS source density with 
reliable shape measurements and photo-zs (Section 4.1) is only of 
order ~45 galaxies per arcmin 2 , thus we cannot directly generate 
predictions for the expected WFIRST source density. Instead, for 
WFIRST, we will use 45 galaxies per arcmin 2 . Our estimates for 
WFIRST will be conservative in this regard. 

As mentioned earlier in Section 5.6, our predictions make sev¬ 
eral simplifying assumptions. First, we use the same shape noise 
as COSMOS (we do not attempt to adjust the measurement error 
component of the shape noise). Secondly, we use photometric red- 
shifts from COSMOS, ignoring the fact for example, that WFIRST 
will detect galaxies in redder bands, and thus a somewhat higher 
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Figure 9. Upper panels: similar to Fig. 7 but here we show only source galaxies that pass also the photo-z cut. Red dots represent all source galaxies with 
both shape and photo-z measurements. Green dots represent only source galaxies that pass our proximity cut with / = 1. Similarly, blue dots represent only 
source galaxies that pass our proximity cut with / = 1.3. Lower panels: Al pa ; rs (< r) divided by A r w ; t hout cut where /V w ithout cut is the extrapolation from 
larger scales of the source counts without any proximity cuts (red points). With f = 1, the overall source density is reduced by 20 per cent. With / = 1.3, the 
overall source density is reduced by 30 per cent (also see Table 2). Proximity effects start to impact :V p;l j rs (< r) at 40 kpc. At the equality radius, the source 
density is reduced by 60-80 per cent for / = 1. 


redshifts. Finally, we also ignore the extra effects of smearing by 
a larger PSF (~ a factor of 1-2 larger compared to COSMOS) on 
Adairs ( r )> the number of source galaxies as a function of transverse 
separation. A more detailed study should account for such effects, 
but our goal here is simply to generate a first order prediction for 
the expected S/N. 

We compute Np^' LID {< r ) for Euclid and use 
A£2f MOS « r ) f° r WFIRST. The expected S/N for one weak 
lensing data point at r < R eq is then computed using Equations 18 
and 19. The results are shown in Fig. 11. Similarly to Fig. 10, 
we observe that the S/N for both Euclid and WFIRST reaches its 
maximum for massive galaxies at lower redshifts. Fig. 11 indi¬ 
cates that Euclid and WFIRST may detect AS on very small radial 
scales with very high S/N greater than 20 for lens galaxies with 
log 10 (M*) > 10.4. Note that this is lower mass than the mass 
probed by strong lensing, which is typically log 10 (M*) > 11.0 
(e.g., Auger et al. 2010b; Oguri et al. 2014; Sonnenfeld et al. 2014). 

Until now, all analyses presented in this section involve the 
photo-z cut, which uses best-fit photo-z on both lens and source 
samples. However, as we observed in Fig. 7, the confusion in source 
selection due to photo-z errors might be buried in our analyses even 
after the close-pair cut is applied (e.g., galaxies that are correlated 


with lens samples). To investigate such possibility, we employed 
several stringent photo-z cuts to separate lens and source galaxies 
and re-compute the predicted S/N. For example, to remove galaxies 
that are correlated with lens samples but have 2 > 2 i, max , we ap¬ 
plied an additional buffer redshift to the maximum photo-z in each 
lens bin as 


z s > Z\ ,max + 0.1(1 + Z\ ,max ) • (20) 


The buffer of 0.1 is motivated by results from Ilbert et al. (2009), 
who report a one-sigma error on the best-fit photo-z of about 
(j z /{ 1 + z) = 0.07 for objects with i + ~ 25.5, which is more 
conservative than the expected errors in Euclid and WFIRST (see 
Laureijs et al. 2011 and Spergel et al. 2013 respectively). The re¬ 
sultant S/N is shown in Fig. 12. The predicted S/N decrease up to 10 
per cent at the higher mass end but the stringent cut (Equation 20) 
does not alter the overall trend - the S/N is maximized at massive 
galaxies at lower redshifts. 
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Figure 10. Predicted S/N for AE measured within R eq as a function of lens stellar mass and redshift in the COSMOS field. Left: predicted S/N for three lens 
redshift bins assuming our fiducial close-pair cut with / = 1. Grey vertical lines represent our lens stellar mass bins. Note that the S/N with z £ [0.3, 0.5] 
decreases rapidly at the highest stellar mass bin because the number of lens-source pairs decrease at high mass end (See Table 3). Right: expected S/N when 
we make more or less conservative cuts on source galaxies in close pair configurations (only the middle redshift bin is shown here). A larger value of f 
corresponds to a more conservative selection of source galaxies. Varying / between 0.8 and 1.3 only has a relatively minor impact on the predicted S/N. 


Euclid WFIRST 




logio( M *[ M 0 D l°g lo (M*[M 0 ]) 


Figure 11. Similar to Fig. 10, but for future surveys - the left panel is for Euclid and the right panel is for WFIRST. Again, we assume our fiducial close-pair 
cut with / = 1 and grey vertical lines indicate our lens stellar mass bins. Short grey thin lines represent the division of one the stellar mass bins into 5 smaller 
bins of width 0.074 dex. Our prediction signal in Figs. 13 and 14 will adopt this finer binning scheme. Note that the S/N with z £ [0.3, 0.5] decreases rapidly 
in the highest stellar mass bin due to smaller number of lens-source pairs (also seen in Fig. 10). The yellow shaded region with log 10 (M*)[MQ] > 11.0 
indicates the typical stellar mass range for strong lensing samples from the SLACS, SL2S, and BELLS surveys with z £ [0.1, 0.9]. Although WFIRST will 
have a higher source density, its smaller area results in smaller S/N compared with Euclid. 
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Euclid 



logio( M *[ M 0 D 


WFIRST 



Figure 12. Similar to Fig. 11 but with a stringent photo-z cut applied (see Equation 20). Compared with Fig. 11, the S/N at the higher mass end decreases up 
to 10 per cent but the overall dependence on the stellar mass and the redshift of lens galaxies does not change. Therefore, the confusion in source selection due 
to galaxies that are correlated with lens samples can be mitigated by such photo-z cuts. This may reduce S/N by a few per cent, but does not alter our main 
conclusion - massive galaxies at lower redshifts yield higher S/N. 


6.3 Predictions for radial bins 

Figs. 10 and 11 show the expected overall S/N for one bin at 
r < R eq . We now consider the expected S/N in finer radial bins. 
Figs. 13 and 14 show the predicted errors on AS for lens galax¬ 
ies with the mean stellar mass of 10 1103 Mq with the bin width 
of 0.074 dex (5 times smaller than our bins in Table 3) at a mean 
redshift of z = 0.68. The radial binning scheme is arbitrary - we 
opt to make 451 bins with logarithmically equal width from 0.013 
to 419 kpc. An NFW profile is assumed for the dark matter compo¬ 
nent and a Hernquist profile is assumed for the stellar component. 
This prediction uses source galaxies that pass all the quality cuts 
introduced in Section 4.1 including our fiducial close pair cut with 

/ = i-o. 

As a reference, in Figs. 13 and 14 we show how AS varies 
when the stellar mass of the lens sample varies by ±0.2 dex (a 
factor of 1.5). As can be seen from these figures, both WFIRST 
and Euclid should be able to tightly constrain a combination of the 
inner dark matter slope and the total stellar mass of galaxies down 
to log(M*) > 10.4. 


7 SUMMARY AND CONCLUSIONS 

Current constraints on dark matter density profiles from weak lens- 
ing are typically limited to radial scales greater than r ~ 50-100 
kpc. Below this scale, there is a paucity of sources (background) 
galaxies due to poor image quality (“seeing”), complicating effects 
such as isophotal blending that inhibits shape measurements, and 
a lack of source galaxies simply due to the relatively small areas 
covered by current lensing surveys. Pushing weak lensing measure¬ 
ments down to smaller radial scales, however, would open up the 
exciting possibility of probing the very inner regions of galaxy/halo 
density profiles. Such measurements are currently only possible for 


strongly lensed galaxies which are rare, occur mainly for massive 
galaxies, and may also represent a biased sample. Weak lensing on 
these small scales would extend constraints to a much wider range 
of redshifts, stellar masses, and to an unbiased selection of galaxy 
types. 

In this paper, we use a weak lensing catalogue from the COS¬ 
MOS survey to investigate the expected S/N of stacked weak lens¬ 
ing measurements on radial scales of only a few tens of kpc. On 
these scales, in addition to dark matter, the weak lensing signal 
is sensitive to the baryonic mass of the host galaxy. Thus, future 
weak lensing measurements will offer the exciting possibility of 
directly measuring mass-to-light ratios and providing independent 
constraints on the Stellar Initial Mass Function (IMF) and on the 
interplay between baryons and dark matter. 

We introduce the “Equality radius”, _R eq , as the radius where 
the dark matter component and the stellar component contribute 
equally to the weak lensing observable, AE. Using a stellar-to-halo 
mass relation that has been calibrated for COSMOS data, we com¬ 
pute the evolution of R eq as a function of lens stellar mass and 
redshift. We show that R eq is of order 7 kpc for M* = 10 9,6 Mg 
and 34 kpc for M, = 10 11,4 Mq at 2 = 0.4. For lens galaxies 
with M» > 10 1 o ' 3 Mq, this equality radius is about a factor of two 
larger than the size of the Kron ellipse which encompasses roughly 
90 per cent of the light. This leaves a very narrow window (width of 
order 10 kpc) with which to measure weak lensing signals on scales 
where the stellar component dominates the lensing signal. The area 
of this window (per lens galaxy) A e g varies between 2 x 10 -4 
arcminute 2 and 0.3 arcminute 2 . We show that A e ff is maximal for 
high mass galaxies at low redshifts due to the fact that a) R eq is 
large for massive galaxies and b) a fixed value of i? eq in physical 
units corresponds to a larger angular size at low redshifts. 

Using COSMOS, we calculate the number of lens-source pairs 
as a function of transverse separation 77 pa i ra (r) down to r = 5 
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Figure 13. Predicted S/N on the A£(r) profile of lens galaxies with log 10 (M*) = 11 for WFIRST. Note that the vertical scale is multiplied by R to 
highlight the error bars. The dark matter profile (AE^m) is shown by the blue line and the stellar component (A£ ste n ar ) is shown by the red line. The black 
line represents the sum of AE^m and AE ste u ar . The green thin vertical line indicates the radius that encompasses 90 per cent of the flux of the lens galaxies. 
The green thick vertical line shows R eq . Grey dashed lines show AE when M* is varied by ±0.2 dex (roughly corresponding to the current uncertainty on 
the IMF). Yellow points indicate radial bins where there are no source galaxies in our COSMOS catalogue. 
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Figure 14. Same as Fig. 13 but for Euclid. Euclid will have a lower source density compared to WFIRST but will observe an area 6 times larger (15000 deg 2 
versus 2500 deg 2 ). This results in a factor of 2 in the expected S/N for Euclid compared to WFIRST. Note, however, that we have neglected other effects (e.g., 
the PSF size of EUCLID is twice as large as WFIRST) which will probably affect these predictions at a similar level. 
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kpc. We show that in the COSMOS catalogue, more galaxies in 
close-pair configurations are rejected because of cuts related to the 
quality of shape measurements rather than due to cuts related to 
the availability of a photo-z measurement. This test is simplistic 
in the sense that we only consider the availability of a shape or a 
photo-z measurement - not the quality of these measurements. This 
question obviously needs to be studied in greater detail, but this 
simple tests suggests that shape measurements rather than photo-z 
measurements may be the limiting factor for these types of studies. 

We investigate how blending and proximity affect the source 
counts as a function of transverse distance from the lens sample. 
We show that within R eq , the number of source galaxies is re¬ 
duced by ~ 20 per cent at AT* = 10 11 o Mq and by 60 per cent 
at M* = 1 O 1 o " 5 M 0 due to blending. This sharp decrease in the 
number of source galaxies on small radial scales is due to the ef¬ 
fects of masking/blending by foreground galaxies (Simet & Man- 
delbaum 2014). We quantify how often blends occurring at R eq 
are due to the lens sample or to other (non lens) galaxies that are 
source galaxies and spatially correlated with the lens sample. At 
R e q at z = 0.67, we show that almost all blends occur with the 
lens sample with M* = 10 10 ' 3 Mq whereas only about 30 per cent 
blends occur with lens galaxies with M* = 10 11 ,o Mq. 

Using a simple criterion based on the overlap between Kron 
ellipses, we study how N pa ,i TS (r) varies with /. Over all the stellar 
mass range, larger / rejects more source galaxies and results in 
lower S/N, which modifies the amplitude of S/N a factor less than 
2 from / = 0.8 to 1.3 and does not change the peak of S/N as a 
function of the stellar mass. 

Finally, we use our COSMOS weak lensing catalogue to make 
a first prediction for the S/N for weak lensing measurements at 
r < Re q for Euclid and WFIRST. Our predictions show that future 
experiments should have enough source galaxies on small scales to 
detect weak lensing down to a few tens of kpc with high S/N. Our 
predictions show that this idea is worth pursuing further - but make 
a number of simplifying assumptions that need to be investigated 
in further detail. 

Our work has focused primarily on quantifying the effects on 
blends on the number of source galaxies at r < R eq . In terms 
of raw numbers, our work shows that even after making a fairly 
conservative selection for our source galaxies, Euclid and WFIRST 
should still have a sufficient number of pairs on small scales to 
measure weak lensing with high S/N. This type of measurements 
however, will face a number of challenges which still need to be 
investigated. The main challenges that remain to be tackled will be 
understanding how to measure both shear and photometric redshifts 
in an unbiased fashion in close-pair configurations. In addition, fur¬ 
ther work is needed to test and calibrate shear measurements in this 
intermediate regime where the shear takes on values above 0.05. 
Also, here we have only considered stellar mass, but in addition the 
gas component may also be non-negligible especially at higher red- 
shifts (Papastergis et al. 2012). Finally, but not least, understanding 
how to separate source galaxies from other galaxies spatially asso¬ 
ciated with lens galaxies and how to correct for boost factors (Man- 
delbaum et al. 2005) represent another challenge to be investigated 
in greater detail. 
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